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Abstract 

We demonstrate the QCD factorization for inclusive hadron production in pA collisions in the 
saturation formalism at one-loop order, with explicit calculation of both real and virtual gluon 
radiation diagrams. The collinear divergences associated with the incoming parton distribution of 
the nucleon and the outgoing fragmentation function of the final state hadron, as well as the rapidity 
divergence with small-x dipole gluon distribution of the nucleus are factorized into the splittings 
of the associated parton distribution and fragmentation functions and the energy evolution of the 
dipole gluon distribution function. The hard coefficient function is evaluated at one-loop order, 
and contains no divergence. 
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1. Introduction. Gluon saturation at small-x (the longitudinal momentum fraction) in 
nucleon and nucleus plays a crucial role in the description of high energy hadronic scat- 
tering [ll-lsj. It has been applied to a wide range of processes which are relevant for the 
current and future experiments at RHIC and LHC [^]. Among these physics, suppression of 
single inclusive hadron production in the forward direction of dAu collisions at RHIC has 
been regarded as one of the most important evidences for the on-set of gluon saturation at 
small-x in a large nucleus (sf. In particular, the dipole gluon distribution applied to hadron 
production in nucleon-nucleus scattering is the same as that in the description of inclusive 
deep inelastic scattering structure function at small-x The duality between these two 
processes, and the universality of the gluon saturation is important to understand the strong 
interaction dynamics in the dense medium [rj. 

However, the experimental data are so far interpreted in the leading order calculation 
in the saturation formalism [8|. In order to solidate the above conclusion, we need to go 
beyond the leading order Born approximation, and evaluate higher order corrections, which 
has not yet been done. A next-to-leading order calculation is highly anticipated in the 
small- a; formalism, to demonstrate the QCD factorization which can be applied to a wide 
range of high energy processes. It has been one of the most important quests in the last 



few years [9|, Il0[. Recent studies on the two-particle production in pA collisions have also 
emphasized the importance of the factorization and universality issues involved in the hard 
processes, in particular, for small-x physics 

In this paper, we will carry out, for the first time, a complete next-to-leading order 
(NLO) calculation for the single inclusive hadron production in pA collisions in the satu- 
ration regime, by employing both the collinear factorization and the high energy small-x 



factorization techniques. Early efforts have been made in Refs. |12l . |13| . Our results will 
not only provide an important estimate of high order corrections, but also pave the way to 
build a QCD factorization for hard processes in the saturation formalism [9]. Moreover, this 



calculation, combining with the associated calculations for inclusive DIS process [10|, will 
help to identify the universality of the dipole gluon distributions in these processes. The 
method developed here will be very useful to other hard processes in small-x physics as well. 
In the process of hadron production in pA collisions, 

p + A^h + X, (1) 

a parton from the nucleon (with momentum p) scatters on the nucleus target (with momen- 
tum Pa)i and fragments into final state hadron with momentum Ph. In the dense medium 
of the large nucleus and at small-x, the multiple interaction become important, and we 
have to perform the relevant resummation to make the reliable theoretical calculations. The 
small-x formalism, the color-glass-condensate or color-dipole formalism, was developed to 
calculate these hard processes We are particularly interested in the scattering process 
with a dilute projectile on a dense target such as the process of ([1]), where the parton from 
the incoming nucleon can be treated as the normal parton distribution and it fragments to 
the final state hadron after multiple interaction with the nucleus target. According to our 
following calculations, the QCD factorization formalism for this process reads as, 

^^^p+A^h+X f dz dx f 

— ^y^2p^ — = z2 J -^—^^f^(^^^^)^h/c{z,fi) J [dx±]Sl^{[x±])na-,c{as,^,[x±]fi) ,{2) 

where ^ = r/xz with r = pxe^ / ^/s, y and p± the rapidity and transverse momentum for the 
final state hadron and s the total center of mass energy square s = {p + PaY, respectively. 
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In the above equation, fa{x) and Dh/c{z) represent the colhnear parton distribution from 
the incoming nucleon and fragmentation function for the final state hadron, where x is the 
momentum fraction of the nucleon carried by the parton a, and z the momentum fraction 
of parton c carried by the final state hadron h, respectively. The response from the nucleus 
target is denoted as S^^{x±) (see the definitions below), depending on the flavor of the 
incoming and outgoing partons and the gluon rapidity Y associated with the nucleus: Y ^ 
ln(l/xg) with Xg being momentum fraction of nucleus entering the hard process. At the 
leading order, they are defined as the two-point functions representing the dipole gluon 
distribution functions in the elementary and adjoint representations for the quark and gluon 
initialed subprocesses Q , respectively. Higher order corrections will have terms that depend 
on the correlation functions beyond the simple two-point functions. Because of this reason, 
the integral [dx±\ represents all possible integrals at the particular order. For example, 
a four-point function (non-linear term of dipole gluon distribution) will contribute to the 
quark channel at one-loop order. The hard factor Tia^c describes the partonic scattering 
amplitude of parton a into a parton c in the dense medium. This hard factor includes 
all order perturbative corrections, and can be calculated order by order. Although there 
is no simple A;^-factorization form beyond leading order formalism [l4|, we will find that 
in the coordinate space, the cross section can be written into a nice factorization form as 
Eq. ([2]). Besides the explicit dependence on the variables shown in Eq. ([2]), there are implicit 
dependences on in the hard coefficients as well. 

The above factorization is derived in the high energy limit, i.e., i/i — i- oo, which is the 
same limit that the dipole formalism was derived for the inclusive DIS structure function 
at small-x [l5|. Two important variables are introduced to separate different factorizations 
for the physics involved in this process: the collinear factorization scale fi and the energy 
evolution rapidity dependence Y. The physics associated with fi follows the normal collinear 
QCD factorization, whereas the rapidity factorization Y takes into account the small-x 
factorization. The evolution respect to fi is controlled by the usual DGLAP evolution. 



whereas that for by the Balitsky-Kovchegov (BK) evolution 16|, Il7|. In particular, 
our one-loop calculations will demonstrate the important contribution from this rapidity 
divergence. 



A different formula has been proposed in Ref . [12[ , where only part of one- loop calculations 
were taken into account, and the rapidity divergence is not identified and the collinear 
evolution effects are not complete. In particular, the formula in Ref. [l2[ was expressed in 
the leading order form, which does not allow higher order corrections. Recently, it has been 
realized that the higher order corrections are important for high p±_ particle production in pA 
collisions, referred as "inelastic" contribution in Ref. [l3|. In this paper, we will demonstrate 
the factorization formula of Eq. (|2]), by a complete one-loop calculations, including not only 



the contributions considered in Refs. [12|, |13|, but also the virtual diagrams which have 



not yet been calculated before. We will take the example of quark channel contribution to 
demonstrate this factorization, and list the results for all other channels in the end. For the 
quark channel contribution: qA q + X, the above factorization formula can be explicitly 
written as 



l^rf^ = j ^-^xq{x,^^)D,/,{z,^^) j {S^ \x^,y^) 



"^(0) (1)- 
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up to one-loop order. The two-point and four-point functions are defined as 



S^'\x^,y^) = —{U{x^)U\y^)) 



N, 



Y 



iV2 



Y 



(4) 
(5) 



IS 



with Nc the number of color in QCD, and U{x_i) = V exp ^^igs /I*^ dx^ T'^A^ 

the Wilson line in the small-x formalism [4] with A~[x~^, x±) being the gluon field solution of 
the classical Yang-Mills equation. For convenience, we also introduce the momentum space 
interpretations of the above functions: J^{k±) = J ^-4j4;T^e~'^''^'^^^^y^^ S'''^\x±,y±), and 



I 



dP'x^cp-y^ifb^ -ik I -{x I -h i )-ik^ , -(b, -y,) 



(2^)4 



S^^\x±,b±,y±). The above two func- 



tions are related as follows: J^{k±) = f (Pki±Q{kj_, ki±). At the leading order in a^, there is 
only two-point function contribution. At the next-to-leading order, a non-linear term of the 
two-point functions will enter. For the gluon channel, a six-point function will apply. The 
goal of this paper is to demonstrate the above factorization formalism, and obtain the hard 
coefficients 7^2^]. In doing so, we calculate the scattering amplitude squared at one- loop or- 
der, and factorize out the divergences associated with the splittings of the quark distribution 
and fragmentation functions as well as the rapidity divergence from the two-point function. 
The hard coefficients are Infra-red and Ultra-violet finite. 

2. One-loop Calculations. The leading order results have been calculated before, from 
which we have 



H 



(0) 
299 



-ik 



-5(1-0 



(6) 



where k± = p±/z and r± = x± — y±. In the following, we will perform the NLO calculations. 
There are virtual and real gluon radiations. We plot the typical diagrams for them in Fig. 1. 
After some algebra, we find that the sum of the virtual diagrams. 



J ^ _ ^ |Cf J d'^qi.1{q±,k_L) + J d'^q±d'^kgi^J{q^,k±,kgi±] 



(7) 



where Cp = {N^ — l)/2Nc, and X and J' are defined as 



X(g_L, k±] 
J{qj_, k±, kgii_] 



q± - k_i 



qi. - iki_ 



[J'{k^)5^^^ {kgl^-k^)-g{k^,kg,^ 



2iq^-^k^)-iq±-k 



gi-L) 



{q± - iki_f{q^ - kgii_) 



In the above results, we notice that the Ultra-violet divergences cancel out among the virtual 
diagrams. This can be seen from the large q± behavior of the above integral and the identity 
relating Q and J-". At large q± limit, the gluonic interaction with nucleus target is not affected 
by the medium effect, and the Ward identity applies, therefore the UV divergence cancels 
out between the self-energy and vertex diagrams. The real gluon radiation diagrams have 



been considered before Il8l | . In order to obtain the single inclusive particle production 
cross sections, we need to integrate out the phase space of the radiated gluon. In the high 
energy limit, this integration will collapse the quadrupole into dipole and non-linear term 
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FIG. 1. Typical Feynman diagrams for the real (a) and virtual (b) gluon radiation contributions 
to the quark production qA q + X at one-loop order. The crosses represent the final observed 
quark, and the blobs in the lower parts of the diagrams represent the multiple interactions with 
the nucleus at small-x. 



of dipoles. At the end of the day, we can cast the real contribution into 



as 
27r2 



+ 



J d'^kg±d'^kgi^J{k^,kg±,kgi±] 



(8) 



8]) when ^ — )■ 1. The rapidity 
(1 - goes 



where x = t/z^. 

There are rapidity divergences in the above results 
divergence appears when the longitudinal momentum of the gluon k^ 
to zero, namely when the rapidity of the gluon goes to — oo. It is important to note that 
the rapidity divergence cancels out for the term proportional to the color-factor Cp- The 
remainder of the rapidity divergence is proportional to color-factor Nc/2, which is exactly 
the same as that in the BK evolution equation [l^, 17 1. It is also interesting to note that the 
rapidity divergence disappears when one integrates over the transverse momentum k± [l9| . 

There are also coUinear divergences from both real and virtual diagrams. We use di- 
mensional regularization (D = 4 — 2e) and follow the MS subtraction. Since the soft gluon 
radiation has been included into the BK evolution of the unintegrated gluon distribution, 
there is no soft divergence. For example, for the virtual diagrams we have 



Air 



CpJ^ikj 



■ 1 kl^ 
— + In ^ 



^(fc^)ln(l-O' 



+ — d kgi±g{k±,kgi±)\n 



{kgl± - ^fcj 
k', 



(9) 



The collinear divergence is represented by 1/e in the above equation. The same divergence 
appears in the real diagram calculations. By summing up both real and virtual contributions, 
we identify the following divergences. 



asNc 
' 27r2 

2n 



de,' 



d%, 



d^ 



r/z 



[S^'\x^,y^)-S^'\x^,b^,y^)] 



(2 



x±,y±) 



(10) 
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where the sphtting kernel is defined as Vqq{^) 



Obviously, Eq. (fTOj) contains 



three divergences: rapidity divergence, and two collinear divergences. The rapidity diver- 
gence f the first term) can be absorbed into the renormalization of the dipole gluon distribu- 
tion |i, Eol, H By doing so, we introduce the Y dependence in the two-point function, 
from which the BK evolution can be understood by identifying dY = dC,'/{l — ^')- The 
collinear divergences symboled by 1/e in dimensional regularization will be absorbed into 
the renormalization of the quark distribution and fragmentation functions, following the 
usual collinear factorization. After subtracting the above divergences, we obtain the hard 
coefficients 



{2Cf - N,) e-'^^-'^ 



n 



(1) 

Aqq 



-ikj_ ■rj_ 



1+e 



(1 + a In (1-0' 



'21 



e « 



X / 



6^'\b^-yA.) / d\\ 



J2 



(11) 



:i2) 



where cq = 2e ""^ with 7^; the Euler constant, and 



'21 





g-i{i-Ofe±-fe± 







bi. ■ {jb^ - rx) 
bl {^b^-r^f 



1 

6^ 



j^^-ik^-b^ 



(13) 



These hard coefficients do not contain any divergence. The calculations for all other partonic 
channels follow the same procedure, and the hard coefficients are calculated up to one-loop 
order. 

The above results clearly demonstrate that the factorization of ([3]) is achieved. The 
collinear divergences associated with the collinear parton distribution and fragmentation 
functions are factorized, and so does the rapidity divergence associated with the two-point 
function of the nucleus. This is a very important step to prove the factorization beyond the 
leading order in perturbation theory. 

Our results also show that we can write down the differential cross section in a factoriza- 
tion form in the coordinate space. The factorization scale dependence in the hard coefficients 
reflect the usual DGLAP evolutions for the quark distribution and fragmentation functions. 
It is interesting to note that similar dependence (associated with r^) has also been found 
in the transverse momentum resummation formalism derived for the Drell-Yan lepton pair 
production in Ref. 
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However, the hard coefficients in our case do not contain double log- 
arithms, and there is no need for the Sudakov resummation for inclusive hadron production 
in pA collisions. 

Great simplification for the above hard coefficients can be found if we take large Nc limit. 
For example, the last term of Eq. ( ITTj) will drop out. This applies to all other quark involved 
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channels as well. Here, we list these hard coefficients in the large Nc limit, 



H. 



(0,1) 
(1) 



H. 



(0,1) 
239 



, H. 



(1) 



H. 



499 



^(1) 



-AtcN, 



2qq - 



Oe"*'^^ "^ In 



A;, 



-47r 



{xl - y±) - yi.)' 



■W(A;j 



(14) 



where yV{ki±,k2±) = e~*'^i-L'(^^~J'^)~*'^2± {j/±-bx)^ g^^^^ j^o^^g chosen /i = co/r± for the fac- 
torization scale to further simplify the above expressions. The complete n dependence can be 
re-stored by using the DGLAP evolution for the relevant parton distributions and fragmen- 
tation functions. With the McLerran- Venugopalan model for the two-point function |3i] , we 
find that the best choice for /i will be in the order of the saturation scale Qs- For the gluon 
channel gA — )■ gX, a factorization similar to Eq. (|3]) was established, where the two-point 
function in the adjoint representation (at large Nc limit as S^'^\x,y)S^'^\y, x)) appears at 
the LO and NLO, and a non-linear term as S^^\x, b, y)S^'^\y, x) enters at NLO. The relevant 
hard coefficients at large Nc limit read as: H. 



(0) 



22 7j(l) 



H2gg, and 



(1) 



655 



^.(y-^b) [1 - ^(1 - 0]^ 1 x±-y_L 



-y± 



(1-0+ e{x^-y^f (6x 



y±y 



_ {b± - yi_) 



S^'\b^-y^) / d\[ 



J2 



(15) 



respectively. The quark loop contribution has been calculated but not included here. 

The above results can be compared to the complete collinear factorization calculation 
for inclusive hadron production at large transverse momentum. This corresponds to the 
situation discussed in Ref. [l3[. Another important issue is the running coupling effects (see 



e.g., Ref. [2l|) at one-loop order, which has not been discussed in this paper. We will carry 
out these studies, together with the phenomenological apphcations to the RHIC and LHC 
experiments, in a separate pubhcation. 

3. Summary. We have demonstrated the QCD factorization for inclusive hadron produc- 
tion in pA collisions in the saturation formalism. The collinear divergences are shown to be 
factorized into the splittings of the parton distribution from the incoming nucleon and the 
fragmentation function for the final state hadron. The rapidity divergence at one-loop order 
is factorized into the BK evolution for the dipole gluon distribution of the nucleus. The hard 
coefficients are calculated up to one-loop order. In principle, using these hard coefficients 
together with the NLO parton distributions and fragmentation functions as well as the NLO 



small-x evolution equation [22|, |23 



one can obtain the complete NLO cross section of the 
inclusive hadron production in pA collisions. 

These results are very important, not only for the phenomenological applications to the 
inclusive hadron production in pA collisions at RHIC and the LHC where active experiments 
are pursued for the study of saturation physics, but also for theoretically promoting the 
rigorous developments toward a complete QCD factorization in small- a; physics j^. Our 
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results will stimulate further applications of the method, in particular, the factorization 
technique, used in this paper to other hard processes involving big nucleus and small- a; gluon 
distributions, as well as those in hot/dense medium. We expect more exciting developments 
along this line. 
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